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Abstract 

We consider atomic chains with nearest neighbour interactions and study periodic and homoclinic 
travelhng waves which are caUed wave trains and solitons, respectively. Our main result is a new 
existence proof which relies on the constrained maximisation of the potential energy and exploits the 
invariance properties of an improvement operator. The approach is restricted to convex interaction 
potentials but refines the standard results as it provides the existence of travelling waves with uni- 
modal and even profile functions. Moreover, we discuss the numerical approximation and complete 
localization of wave trains, and show that wave trains converge to solitons when the periodicity length 
tends to infinity. 
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1 Introduction 



Wc consider infinite chains of identical atoms with unit mass that are coupled by nearest neighbour 
interactions. The dynamics in such chains is governed by Newton's equations 

where Xj (t) denotes the position of the jth atom at time t, and (f> is the interaction potential. Restating 
(1) in terms of atomic distances rj{t) = — Xj{t) and atomic velocities Vj{t) = Xj{t) we find 

r,{t) - v,+,{t) - v,{t), v,{t) - <^>'{r,{t)) - <J''(r,-i(t)). (2) 
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In this article we allow for arbitrary convex interaction potentials $ and refer to (1) as Fermi-Pasta-Ulam 
(FPU) chain although the potential in the original paper [FPU55] was a quartic polynomial. 

FPU chains can be viewed as simple toy models for crystals and solids and allow to study some 
essential properties of nonlinear clastic materials. Even though (1) is a strong simplification of a real 
material it obeys a very complex behaviour and currently we are far from a complete understanding of 
its dynamical properties. 

During the last decades, a lot of research addressed the existence and properties of travelling waves in 
FPU chains because they can be viewed as elementary waves and provide a lot of insight into the energy 
transport in nonlinear media. Such travelling waves solve nonlinear advance-delay differential equations 
and it is a fundamental mathematical problem to characterise the solution set to those equations. 

A further motivation for the study of travelling waves is related to atomistic Riemann problems and 
self-thermalisation of FPU chains: Starting with piecewise constant initial data for the atomic distances 
and velocities, solutions to (2) are self-similar on a macroscopic scale and involve dispersive shocks, that 
are fan-like structures with strong microscopic oscillations, see Figure 1. It is known from the theory 
of integrable systems and numerical simulations that the oscillations within a dispersive shock can be 
described by modulated travelling waves, compare [FV99, DII08] and references therein. 
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Figure 1: The atoms in a dispersive shock self-organize into modulated travelling waves. Left. Snapshots 
of atomic distances and velocities against the scaled particle index for different macroscopic times. Right. 
Corresponding family of wave trains with super-sonic soliton as 'homoclinic limit'; the picture shows 
the density plots of six local distribution functions in the (r, u)-planc. Modulation theory predicts that 
each of these local distribution functions corresponds to a travelling wave whose parameters depend on 
macroscopic time and particle index. 




Travelling waves and main result Travelling waves are exact solutions to (2) and satisfy the ansatz^ 

r,{t)^R{^), v,{t)^V{'f), ^^kj+iot, (3) 

with phase (p, wave number k, (negative) frequency to, and profile functions R and V. Inserting (3) into 
(2) we obtain the nonlinear advance-delay differential equations 

^4^R-ny, ^^v^v^nR), (4) 

where and denote the forward and backward difference operators with shift fc, respectively. 
Depending on the properties of the profile functions we distinguish the following cases: 



^The ansatz (3) is slightly more general than the usual one, which assumes that the atomic positions Xj depend on the 
phase variable (p. In fact, the ansatz Xj(t) = X{kj + Lot) is not invariant under Galilean transformations x x + v^t, 
whereas (3) has this property as it respects r ~^ r,v v + vq. 
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1. Wave trains or periodic waves: R and V are periodic, 

2. Solitons or homoclinic waves: R and V arc localized over a constant background state, 

3. Fronts or heteroclinic waves: R and connect different constant background states, 

4. Oscillatory fronts: R and connect different asymptotic wave trains. 

Notice that our usage of 'soliton' is quite sloppy: Localized travelling waves are sometimes called 'soli- 
tary waves', and 'soliton' then refers to a solitary waves that survives collisions with other such waves 
unchanged. 

In this article we show the existence of wave trains and solitons with unimodal and even profile 
functions R and V , where even means as usual invariance under Lp —1^9, and unimodal functions are 
monotone for both (p < and (p > 0. Our main result can be stated as follows. 

Theorem 1.1. Under natural regularity assumptions on the convex potential $ there exists a four- 
parameter family of wave trains, and if $ satisfies additionally some super- quadratic grow conditions 
then there exists also a three-parameter family of solitons. Moreover, the profile functions R and V for 
both families are unimodal and even. 

Close to our work are [FW94, SW97, FV99, PPOO, SZ07] where the existence of travelling waves is 
likewise studied in a variational framework. Therefore we shall compare both our method (§2.2) and 
results (§4.3) with those presented therein. We also refer to the numerical study [DEFW93], to [FM03] 
for existence results in 2D lattices, and to [looOO] which proves the existence of small-amplitude travelling 
waves by means of center-manifold reduction. Moreover, the existence of fronts is studied in [HR08a], 
and [SZ08] concerns oscillatory fronts in FPU chains with bi-harmonic potentials. 

All results presented below concern solely wave trains and solitons in FPU chains, but the method 
can also be applied to other Hamiltonian lattices with convex potential energy V as for instance Klein- 
Gordon chains with convex on-site potential, see [Her09] , and atomic chains with next-nearest neighbour 
interactions. 

We emphasize that we are not able to provide uniqueness results for travelling waves. Uniqueness of 
relative equilibria in Hamiltonian lattices is a notoriously difficult problem, and almost nothing is known 
about it. The only available results concern either the near-sonic limit [FP99], or systems where the 
travelling wave equation can be solved explicitly. Examples are the Toda chain [TodSl], the discrete 
nonlinear Schrodinger equation [HLM94], the harmonic chain, and the hard-sphere model [DHM06]. 

Overview on the proof and organisation of the paper In a first step we reformulate the travelling 
wave equation (4) in terms of a normalised profile function W E L^. More precisely, we show that (4) 
can be transformed into a nonlinear eigenvalue equation 

uj^W = dViW), (5) 

where V(W) is the potential energy of a travelling wave. The profile function W has no physical meaning, 

but determines R and V via Rif ) = + J^'^ ' W{(f) d^p and V{(p) = vq + LjW{ip), where rg and vq are 
suitable normalization constants. 

Our approach relies on a combination of variational and dynamical concepts, and can be summarized 
as follows. 

1. Equation (5) is the Euler-Lagrange equation for the optimisation problem V{W) — *■ max subjected 
to the constraint W G Bj, where 7 is a free parameter, B.y C denotes the ball of radius ^/2^ , 
and is the Lagrangian multiplier. 

2. There exists an improvement dynamics W t-^ T^lW] on Bj that increases the potential energy. 
Moreover, each stationary point of this dynamics solves (5), and vise versa. 
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3. There exist non-trivial cones S that are invariant under the improvement dynamics. Consequently, 
each maximiser oi V in S H B-y is a travelling wave, see Theorem 2.3. 

We emphasize that the convexity of $ is essential for this approach as it is intimately related to both the 
properties of T-y and the existence of non-trivial invariant cones. 

A major part of the mathematical analysis done in this paper is needed to show that there exists 
a maximizer of V in S-y. In the wave-train setting we can use rather simple compactness arguments as 
the functional V is continuous with respect to the weak topology in L^. In the soliton setting, however, 
we lack the weak compactness of 7^, and the existence proof for maximizers requires more sophisticated 
arguments. Our main technical result in this context is Lemma 4.7, which implies that (for certain S) 
the maximising sequences for V in are localized, and hence pre-compact in the strong topology. 

This paper is organised as follows. In §2.1 we derive the fixed point equation (5) for both wave 
trains and solitons. The details of our variational approach are presented in §2.2, and §2.3 concerns the 
analytical properties of the underlying functionals and operators. In §3.1 wc continue with the existence 
proof for wave trains and present some numerical simulations; the complete localisation of wave trains 
is studied in §3.2. The existence proof for solitons is contained in §4.1 and relies on a natural condition 
for the super-quadratic growth of the functional V. In §4.2 and §4.3 we then discuss the corresponding 
properties for the interaction potential $. Finally, inspired by the notion of F-convergence we show in 
§4.4 that wave trains converge to solitons when the periodicity length tends to infinity. 

2 Variational approach 

In this section we transform the travelling wave equation into a fixed point equation for a normalised 
profile function W, and describe our variational approach to existence results for both wave trains and 
solitons. To point our the key ideas we start with more formal consideration in §2.1 and §2.2, and 
postpone the analytical details to §2.3. 

2.1 Travelling waves as eigenfunctions of nonlinear integral equations 

In what follows we assume that the periodicity length of wave trains is given by 2L with < L < oo, 
and regard the corresponding profile functions R and V as defined on [—L, L]. Moreover, we identify the 
soliton case with L = cxj by considering R and V as functions on [— oo, oo], or, equivalently, as functions on 
the Alexandrov compactification of R. In other words, in both cases we impose the 'boundary conditions' 
R{L) = R{-L) and V{L) = V{-L). 

In what follows we denote by 1?{[~L, L\) with L e (0, cx)] the Lebesgue space of all square- integrable 
functions on [— L, L]. and if there is no risk of confusion wc write instead of L^([— i, L]). 

Our first aim is to transform the travelling wave equations for wave trains and solitons into eigenvalue 
equations for certain nonlinear integral operators defined on L^. For this purpose we define two linear 
averaging operators, and normalise the potential $. More precisely, for given reference distance tq we 
define the potential ^,-0 by 



which is normalised via ^'ro(O) = ^ro(O) ~ ^-^^id $"-,(0) = $"(ro), and mention that this normalisation 
respects the convexity of $. The eigenvalue equations for solitons and wave trains involve the averaging 
operators 



ip+k/2 





ip-k/2 
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and 

L L 

AkW = AkW J AkW{>p) dip = AkW - ^ y W{>p) d>p. 

Notice that the operator Ak is well defined and symmetric on for both finite and infinite L, compare 
Lemma 2.5 below, whereas Ak is well defined for L < oo only. 

Wave trains and normalisation via mean values In order to reformulate the travelling wave 
equation (4) for L < oo we introduce the mean values of wave trains by 



1^ ±j 



L 

_ 1 

r.. - 

-L 

Observation 2.1. With the identification 

i?(^-fc/2) = rav + (AVK)(^), Viip)=v,,+LuWiip), (6) 
for some profile function W with J^j^ ^if) 'if ~ the integral equation 

uj^W = AkKJ^kW). (7) 
is equivalent to the wave train equation (4). 

Proof. First suppose that R and V solve (4), and let W ~ oj^^{V — Wav)- Then, the first equation in (4) 
implies 

-^Riv - k/2) = W{ip + fc/2) - W{p - k/2) = (AkW^ (if), 

and hence R{p — k/2) ~ {AkW){(p) + C\ for some constant ci. Integrating this with respect to p wc find 
Tw = + ci; and hence (6)i. Moreover, the second equation in (4) provides 

w'^W^M = *'(rav + {AkW){-p + fc/2)) - $'(rav + {AkW){p - fc/2)) 

= <i>:.^^ ({AkW){p + fc/2)) - <i>:,^^ ({AkW){p - fc/2)) , 

and integration with respect to p) gives uj^W = Ak^'r._^^{AkW) + C2 for some constant C2. The condition 

■jj- J^j^ W{p) dp ~ now implies C2 = and hence (6)2. Now let a solution W of (7) be given, and define 
both R and V by (6). Then, the first equation from (4) holds by construction, and the second one follows 
from (7) by differentiation with respect to p. □ 

Solitons and normalisation via background states For solitons it is much more convenient to base 
the normalisation on the background states 

rbg := lim R{p), Vhg := lim V{p), 

if — >±OC — >±00 

because then the normalised profiles are localised. 
Observation 2.2. With the identification 

R{p-k/2)=ry,g + {AkW){p), V{p)=v^g + LoWip), (8) 
for some profile function W with lim^^±oo W{p) = the integral equation 

iu^W = AkK^^iAkW). (9) 

is equivalent to the soliton equation (4). 
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Proof. Similar to that of Observation 2.1. 



□ 



We can consider the integral equation (9) also for L < oo, and readily verify that each solution defines 
a wave train via (8) . This normalisation for wave trains seems to be artificial because for finite L we have 
no immediate interpretation of the parameter rbg. However, in §4.4 we rely on this setting and show that 
the L-periodic solutions Wl to (9) converge to solitons as L ^ oo, see Corollary 4.21, an this implies 
Rl{±L) ^ rbg. 

Parameter dependence of travelling waves The travelling wave equation (4) obeys a simple scaling 
symmetry. In fact, the solution set to (4) is invariant under 

R{ip)-^Ri\(p), V{ip)-^V{X(p), k^X-^k, Lu-^X-'^uj, L-^X-^L, (10) 

with arbitrary A > 0, corresponding to 

W{ip) --^ XW{Xip) with rbg~^rbg, Wbg~~^Wbg, rav~^rav, Wav-^Wav 

Up to this scaling there remain four independent parameters for wave trains. A natural choice for 
nonlinear potentials is to fix the length parameter L, and to regard Tav, i^av, k, and uj as independent 
parameters, see [DHM06]. In this paper, however, we prefer to fix the wave number k, so that wave trains 
are parametrised by rav, Vav, L, and the phase speed a := uj/k. Accordingly, we parametrize solitons by 
'^bg, Wbg and a. 

Notice that the velocity parameters v^-v and Ubg do not appear in the travelling wave equations 
due to the Galilean invariance of FPU chains. Both Uav and Ubg are thus trivial parameters, which 
however become important when studying modulated travelling waves for which the parameters vary on 
a macroscopic scale, compare [DIIM06, DH08]. 



2.2 Variational structure for wave trains and solitons 

In view of the scaling and reformulation results from §2.1 we start to simplify our setting. In what follows 
we consider an arbitrary convex and smooth potential <& normalised by <I'(0) = $'(0) = 0. Moreover, we 
restrict to = 1, so that the frequency lo equals the phase speed cr, and consider solely the following two 
averaging operators 

¥>+l/2 L 

{AW){^)= I W{^)d^, {AW){^) = {AW){^)^^ I W{^)d^. (11) 



Moreover, in order to point out the similarities between wave trains and solitons we refer to an abstract 
averaging operator A, which equals either A or A, and consider the general travelling wave equation 

a^W ^ A<P'iAW), (12) 

with W G L^([— L, L]), where L might be finite or infinite. In what follows we call W a travelling wave 
if and only if there exist cr^ > such that (12) is satisfied. 



Starting point for the variational formulation is the observation that (12) is the Eulcr-Lagrange equa- 
tion to the action functional 

C{a^, W) = /C(o-2, W) - V{W), 

where 

L L 

IC{a^,W) = ^a'^ J Wiipfdip and P = J <i>{{AW){ip)) dip (13) 
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are the kinetic and potential energy, respectively. 

The first rigorous existence result for wave trains and solitons under super-quadratic growth assump- 
tions for $ was given by Friesecke and Wattis [FW94]. The key idea is (in our notations) to minimise 
the L^-norm of W under the constraint of prescribed potential energy, where l/cr^ plays the role of a La- 
grangian multiplier. The existence of corresponding minimisers was then established by means of Lions' 
concentration-compactness principle [Lio84] . 

Smets and Willem [SW97] prove the existence of solitons by showing that for super-quadratic $ the 
functional £ satisfies the assumptions of (a modified) Mountain Pass Theorem. Recently, these results 
could be improved by Schwetlick and Zimmer. In [SZ07] they require the super-quadratic growth to hold 
only asymptotically, so certain double-well potentials are admissible. 

Similarly to [SW97], Pankov and Pfliiger [PPOO] apply the Mountain Pass Theorem to wave trains, 
and pass to the limit L ^ oo by means of concentration compactness. Moreover, they present a different 
existence proof for solitons which is based on the Nehari manifold of C. 

A new idea for convex potentials $ was introduced by Filip and Venakides [FV99] by proposing to 
maximise the potential energy V under the convex constraint W £ Bj with 

B-, = {WeL^ : i||T^||B7}- 

A first advantage of this approach is that L < oo implies the functional V to be continuous with respect 
to the weak topology in L^, so that the existence of wave trains follows from elementary principles of 
infinite-dimensional convex analysis. Secondly, the improvement operator 7^ appears naturally in this 
context, and gives rise to effective approximation schemes for wave trains, see [Her05, DH08] and the 
numerical simulations below. 

Our method is also based on the constrained maximization of the potential energy but yields improved 
results as it exploits several invariance properties 7^. 

Constrained maximisation of the potential energy via the improvement operator 

In view of (12) we formally define the improvement operator by 

where the operator dV is the Gateaux derivative of V, that means dVlW] = A^'{AW) for all W e L^. 
By construction, each fixed point of 7^ is a travelling wave with ^ || \\2— 7, and vice versa, where 
the speed is given by ct^ =|| ||2 /V27. 

In the next Section §2.3 we exploit the convexity of $ and derive the following building blocks for the 
existence proof. 

1. dV respects the positive cone 

U :={]¥ eL^ : W{~ipi) = W{ipi) and W{ipi) > W{ip2) for almost all < ipi < ip2 < L}, 

which consists of all functions on [— L, — L] that are even and unimodal. Moreover, for A = A the 
operator dV also respects 

Af ■.= {W el^ : W{ip) > for almost all if e [-L, -L]}, 

which is the cone of all non- negative functions. 

2. 7^ is well defined on \ A4 and maps into dB-^ \ A4, where 

M:^{W eL^ : V{W) = O} 
is the set of all global minimisers of V. 
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3. increases the potential energy, that means V{TyW) > V{W) for all W ^ M, where equahty 
holds if and only if = Ty[W]. 

Wc arc now able to describe the key principle that provides the existence of travelling waves. 

Theorem 2.3. Let 5 C 6e some positive cone that is invariant under the action of the operator dV . 
Then, the set S^\M. with = S C]B^ is invariant under the action oj T^,and each proper maximiser of 
V in is a fixed point of Tj, and hence a travelling wave. 

Proof. The invariance of is implied by the assumption on S and the properties of 7^. Now let W be 
a proper maximiser. Then, riTy[W]) = V{W) > minP|5^ implies both Ty[W] ^ W and W ^ M , a.iid 
we conclude that W is in fact a travelling wave with cr^ > 0. □ 

In what follows the cone S is given by either U oi U M. Since these cones are not open in V? the 
fact that each maximiser must satisfy the Euler-Lagrange equation (12) with multiplier cr^ is not clear a 
priori but provided by the invariance of under 7^. 

Theorem 2.3 yields only a sufficient condition for the existence of travelling waves. In fact, to show 
that V attains its maximum in S-y is not trivial at all (at least in the soliton case), and requires a better 
understanding of the energy landscape in S-y. In our analysis we follow the direct approach and show 
that maximising sequences for V are compact in some appropriate topology in L^. More precisely, for 
wave trains we use weak compactness, whereas in the soliton setting we establish the strong compactness 
for maximising sequences. 

2.3 Some pieces of functional analysis 

Here wc prove the aforementioned properties of the improvement operator 7^. To this end wc rely on 
following standing assumptions on the potential $. 

Assumption 2.4. For given j > we assume that the interaction potential $ has the following properties 
on the interval [— y^27, \/2^\. 

1. (Smoothness) $ is at least C^, 

2. (Convexity) $" > 0, 

3. (Normalisation) = $(0) = $'(0) and $"(0) /3 > 0, 

4. (Non-triviality) $ does not vanish identically. 

The restriction to the interval [—^/2^, is natural in our context, because W G implies 

1 1 1 1 00^ •\/27, see Lemma 2.5. As a consequence of Assumption 2.4 we find 



for all r with \r\ < \f2Fi. Moreover, the non-triviality condition implies J\A = {W : 'P{W) — 0} and 
Br^\M. ^ {0}, and hence each maximiser of V in is proper. 

Within this section, the parameter L can take arbitrary values in (0, 00], and and W^^^* with 
1 < p < (X) denote the usual Lebesgue and Sobolev spaces on [— i], where 



< $(r) < ir2(/3 + o(|r|)). 



*'(-k|)<0<$'(|r|) 



(14) 




gives the inner product in L^. 
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Properties of the averaging operators A and A 

We summarize some elementary properties of tlie averaging operators that are used in tlie proofs below. 
Lemma 2.5. For any L, the operator A is well defined on and has the following properties: 

1. A maps into n L°° with \\AW\\oo<\\W\\2 and \\AW\\2<\\W\\2. 

2. A maps into \N^'^ with {AW)' {ip) = W{ip + l/2) - W{ip - 1/2). 

3. A is self-adjoint on L^. 

4- If a sequence (W„)^ converges weakly in to some limit Woo , then {AWn\^ converges strongly in 
L^([— i, L]) for each finite L < oo with L < L. In particular, for L < oo the image of each bounded 
set in under the operator A is pre-compact in with respect to the strong topology. 

5. In the wave train case (L < oo) the operator A is compact. Moreover, the mth eigenvalue 
(m=0, 1, 2, ...) is given by 

£'m = 0(^^), Q{g) := sin{g), 

and the corresponding eigenspace is spanned by cos {^^■) and sin {^^■) 

6. In the soliton case (L ^ oo) the operator A is not compact anymore, because it has continuous 
spectrum speC|_2yl = {©(£') ■ Q G M}. In particular, spcC|_2^^ = [0, 1]. 

Moreover, for L < oo we have A : ^ n L°° with MW^||oo<|| VF||2 and \\AW\\2<\\W\\2. 
Proof. Definition (11) gives 

L 

AV{^) = J xiv- -'^Wis) ds = x{^^-)*V (15) 

-L 

where x abbreviates the indicator function of the interval [—1/2, 1/2], and * denotes the convolution 
operator. Holder's inequality provides |^1^(</3)|^ < ^"^^1^2 ^(*)^ ^^'^ from this we readily derive the 
first assertion. The proofs of the second and the third claim are then straight forward. Now suppose that 
Wn Woo weakly in 1?{[—L, L\), so that AWn AWoo point-wise thanks to (15), and this impHes the 
strong L^-convergence on each finite interval [—L, L] due to the uniform L°°-bounds. 

Towards the spectral properties of A we study how acts on plane waves. A direct calculation shows 
that each plane wave Ek{(p) ^ f ^ e^^'^ satisfies the eigenvalue equation 

AEk = Q{k/2)Ek 

pointwise, and this implies the fourth and the fifth assertion. Finally, for L < oo we have 

AV{ip)=x{ip~-)*V, xi^)=xiv)^^. 

This implies |^T/F(v3)p <|| X l|oo|| xCv' ~ O^lli^ l-^M^(v?)P and in turn the desired properties of A. □ 
As a consequence of Definition (11) and Lemma 2.5 we easily find 

1/2 

keTL2A=\weL^:W{-) = W{- + l), J M^((^) dv? = o|. (16) 

-1/2 

In particular, the kernel of A is trivial if either L is irrational, or L = oo. Moreover, for L < oo we have 

kerL2^= kerL2^©span{l}. (17) 



9 



Properties of the potential energy functional V 

We rely on Assumption 2.4 and use standard methods from convex analysis to prove some properties of 
V and its derivative. All results are formulated in terms of the abstract averaging operator A and hold 
both in the wave train and the soliton case. 

Lemma 2.6. The functional V is well- defined, bounded, continuous and Gateaux-differentiable onB^, 
and its derivative dV ^ Ao 9$ o A is a monotone operator and maps continuously into . Moreover, 
for arbitrary Wi , W2 G B-y we have 

PiW2)~r{Wi) > f \\AW2-AWi\\l +{dV[Wi], W2-W1). (18) 

and 

{dr[W2] - dr[Wi], W2 ~Wi)>m \\AW2{ip) - AWi{ip) \\l (19) 
where the monotonicity constant m is given by m = inf $"(r) > f3. 

Proof. For all W in Bj we have ||.4VF||oo< a/27, and Assumption 2.4 implies 

{r)\ < C \r\ , $(r) < iCr^ 
for all r with \r\ < y/2rf, where C = sup|^|<y2^ |<f>"(r)| > (3. Consequently, wc find 

< V{W) < \G II AW ||^< 7 C, II y^9$[y^iy] ||p< CP II W ||p< CP |1 W ||p, 

and all assertions concerning the continuity and boundcdncss of both V and dV follow immediately. 
Now let W\, W2 S B.y be fixed, and notice that the convexity inequality (<E>'(r2) — ^'{ri)){r2 — ri) > 
m{r2 — ri)'^ with = Wi{Lp) implies (19) by integration w.r.t. ip. To prove (18), let 77 £ [0, 1] and 
consider W {rj) := (l — rj)Wi + riW2 G By as well as 

p{v) ■.= V{W{v)). 

The function p is well-defined and diffcrcntiablc with respect to ry, and using (19) wc find 

-^p{^) = {dV[W{r,)], W2-W1) 
dry 

> ir^{dV[W{r])] - dr[Wi], 7JW2 - 7]Wi) + {dV[Wi], W2 - Wi) 
= ri-^{dr[W{ri)] - dr[Wi], W{t]) - Wi) + {dr[Wi], W2 - Wi) 

> T]-^m II AW{t]) - AWi \\l +{dV[Wi], W2 - Wi) 

> rjm \\AW2 - AWi \\l +{&P[Wx], W2 - Wi). 

Finally, wc integrate the last estimate from ry = to r/ = 1, and this gives (18). □ 

The convexity of V implies that each trivial travelling wave with ct^ = must belong to M , the set 
of all minimiscrs of V . 

Remark 2.7. We have dV[W] ^ for all W e By\ M. 

Proof Assume there exists some W € B-f\M with dr[W] = 0. Then, (19) with W2 = and Wi = W 
provides AW = and hence 'P{W) = 0, which is a contradiction. □ 

Notice that for non-degenerate potentials $ with <i>(r) > for all r 7^ we have A4 = ker^ , where 
ker^ is given in (16) and (17). 
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Properties of the improvement operator 

First we show that the cones U and N are invariant under the action of both A and Here again the 
convexity of $ enters as it guarantees that $' increases monotonically. 

Lemma 2.8. The cones lA and J\f are convex, closed under weak and strong convergence in L^, and 
are invariant under the action of dV and A. Moreover, U is also invariant under the action of A (for 
L < oo). 

Proof. The only non-trivial assertion is the invariance of U under A. To prove this, we fix G and 
consider Y = -i^AW with Y{ip) = W{ip + 1/2) - W{ip - 1/2) thanks to Lemma 2.5. This hmction is 
odd as W{lp) — W{~(p) implies 

Y{~ip) = W{-ip + 1/2) - W{~ip - 1/2) = W{ip - 1/2) - W{lp + 1/2) = -Y{ip). 

Hence it remains to show that Y{lp) < for all <p > 0, which is equivalent to 

W{ip + 1/2) <W{ip~ 1/2), ip>0. 

For 1/2 < (y9 < L - 1/2 this estimate follows from < ip - 1/2 < + 1/2. Moreover, for < ip < 1/2 it 
holds thanks to W{-ip +1/2) = W{ip - 1/2) and < -(^ + 1/2 < + 1/2, and for L - 1/2 < 95 < L it 
is a consequence of + 1/2) = Vl^(2L - - 1/2) and 2L - - 1/2 > </? - 1/2. □ 

Lemma 2.9. The operator 7^ maps B-y\ A4 continuously into dB-y \ Ai and satisfies 

r{T^[w])~v{w) > \m \\Ar^\w\~Awfi (20) 

for allW ^B^\M.. Moreover, the equality sign holds if and only if W is a fixed point of 7^ . 

Proof. Lemma 2.6 and Remark 2.7 imply that 7^ is well defined and continuous on Ai. Moreover, 
II Ty[W2] II 2= holds by definition. Now let Wi £ B-y \ M he fixed, and set W2 := T^[Wi] and 

cr| :=|| dr[Wi] II2 /\/27 > 0. Hence, alW2 = dV[Wi], and from (18) we infer that 

r{W2) - V{Wi) - im |j^M^2 - AWi \\l > <^2^{W2, W2 - Wi) 

>^2m\W2\\l - ||Vl^2||2||M/l||2), 

which gives (20) due to || Wi ||2<|| W2 \\2= \f2r^- Moreover, we find an equality sign in the second estimate 
if and only if || W2 \\\=\\ Wi \\l= {Wi, W2) , that means if and only if Wi = W2. □ 

With Lemma 2.9 we have derived all ingredients that we had used in the proof of Theorem 2.3. 

3 Wave Trains 

As a first application of Theorem 2.3 we establish the existence of wave trains in §3.1 and proceed 
with some comments on the numerical computation of wave trains. Afterwards we study the complete 
localisation of wave trains in §3.2. 

3.1 Existence results 

Our first existence result concerns wave trains that are rcnormalised via their mean values. This corre- 
sponds to i < 00, A = A, and S = IA. 

Theorem 3.1. For each L < 00 and 7 > there exist a unimodal and even wave train W such that 
\ \\W\\l= 7 and a'^W = M' {AW) for some cr^ > 0. 

Proof. With respect to the weak topology in , the functional V is continuous and the set S-y =Ur\B-f is 
compact. Hence there exist a maximiser W, which is moreover proper, that means supT'l^^ ~ V{W) > 
= minT'l^^, because $ is non-trivial. The desired result now follows from Theorem 2.3. □ 
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Ill view of Observation 2.1 and the scaling (10) we infer that Theorem 3.1 imphes the existence of a 
four-parameter family of solutions {R, V) to the original travelling wave equation (4) with fixed L. This 
family is parametrised by rav, t'av, fc, and 7, and for nonlinear potentials we can moreover expect (at 
least locally) that 7 can be replaced by oj. 

Similar existence result for wave trains in convex FPU chains arc proven in [FV99, DHM06], but 
provide only W S dB^. Our method improves these results as it establishes the existence of wave trains 
with the additional property W £U, which in turn implies R £U and aV G U. This sheds light on some 
observations from [DH08, HR08b]: The traces of travelling waves found in the numerical simulations of 
initial value problems for (1) typically encircle convex sets in the (r, w)-plane, compare the right picture in 
Figure 1. In particular, these curves have exactly two extrema in both the r-direction and the w-direction, 
and hence they correspond to unimodal profile functions R and V . 

We proceed with some remarks concerning the uniqueness of wave trains. The norm constraint 
\ II II 2= 7 alone is not sufficient for uniqueness as the travelling wave equation is invariant under shifts 
in Lp. Moreover, since the set of all 2i-pcriodic functions, with mL = L for some m S N, is invariant 
under the action of 7^, we can construct a whole family of wave trains satisfying the norm constraint. 

From these considerations we conclude that any uniqueness result for wave trains must prescribe 
further properties of the profile function W . Motivated by numerical simulations we conjecture, that for 
each 7 there exists exactly one travelling wave with W £ U and ^ II 111= 7; but we are not able to 
prove this conjecture. 

Finally, we use the same arguments as in the proof of Theorem 3.1 to derive a similar existence result 
for wave trains in the setting L < 00, A ~ A, and S = U Af. 

Lemma 3.2. For each L < 00 and 7 > there exist a wave train W G U f] M C] dB-y such that 
a^W = A'i''{AW) for some > 0. 

Numerical computation of wave trains It is natural to use the improvement dynamics 

W € ^ Xy[W] € (21) 

for the approximation of wave trains, and a corresponding discrete scheme is readily derived and imple- 
mented. It was proven in [Her05] that the orbits generated by (21) are compact in the strong topology, 
but from a theoretical point of view this result remains unsatisfactory due to the lack of uniqueness. So 
it is neither clear that maximisers of V in are unique, nor that all fixed points of Xy are (global) 
maximisers. 

In numerical simulations, however, we found (21) to have good properties. For a wide class of potentials 
we observed rapid convergence to a unique limit independent of the chosen initial data. In Figure 2 we 
present the numerically computed profiles W for different values of 7 with A = A and 

$(0 = cosh(r) - 1. (22) 

For small 7 we can approximate the potential by the harmonic one $harm(^) = $"(0)r^, and hence the 
profile W is close to a rescaled plane wave. For increasing 7, however, the nonlinearity dominates and 
the profile function becomes tighter. Figure 3 shows the corresponding traces in the (r, 7j)-plane, these 
are the curves 

compare Observation (2.2). Surprisingly we find a nested family of curves, that mean the traces for 
different values of 7 do not intersect but fill out a convex region. We are not able to prove this observation 
but mention that a similar phenomenon occurs when FPU chains generate dispersive shocks, see [DHR06, 
HR08b]. 
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Figure 2: Profile functions W for several values of 7 with L = 2 and $ as in (22). 
TW-traces TW-traces, zoom in 





+0.00 



+4.39 



+0.00 



Distance 



Distance 



Figure 3: Traces for the wave trains from Figure 2. 



3.2 Complete localisation of wave trains 

It is well known for strongly nonlinear potentials that in certain limits the wave trains (and even the 
solitons) localize completely, in the sense that - under a suitable rescaling - the profile functions W 
converge to the indicator function of an interval plus a constant background state. Such profile functions 
are, up to renormalisation. equal to the profile functions of travelling waves in the hard-sphere model for 
the atomic chain, in which all atomic interactions are described by elastic collisions. Thus the effect of 
localisation can often be linked in a natural way to the high energy limit of travelling waves. For more 
details we refer to [Tod81, FM02, Her05, DHM06]. 

In this section we discuss the localisation phenomenon in our context, and aim to derive a localisation 
criterion for wave trains. To keep the presentation simple we solely consider non-negative and unimodal 
profile functions, that means we investigate the localisation of solutions to (12) with L < 00, A ^ A, and 
S = 14 OJV. Moreover, for our purpose it is sufficient to assume that the localised limit profile is given by 



Wch{^) = Xr 1 1,('^) = 

I 2' 2' 

see Figure 4. It is easy to check that this profile satisfies 



if 1^1 
if 1^1 



(23) 



II WcL ||2= 1, (AWclKv) - max{l - |<^| , 0}, P(I^cl) = 2 / $(s) ds 



(24) 



In what follows we consider sequences (<&«)„ of rescaled potentials, where each $„ : [0, 1] ^ M satisfies 
Assumption 2.4. Moreover, we refer to a sequence of profile functions Wn C as a corresponding sequence 
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of maximisers, if Wn is a maximiser of Vn in S1/2 for each n, where Vn is the potential energy functional 
(13) corresponding to <&„. 



WcL{ip) versus tp 



■AWchiip) versus ip 




-1/2 1/2 



-1 -1/2 1/2 



Figure 4: The functions VFcl and ^H/cl- 

We say, a sequence of such potentials ^^^^ complete localisation property on [0, 1] if any 

corresponding sequence of maximisers converges strongly in to Wcl- Our main result in this section 
is a necessary condition for the complete localisation of wave trains and is implied by the following 
observation. 

Lemma 3.3. We have MW^||oo<MM/cl ||oo= 1 for any W G S1/2 with W ^ Wcl- 
Proof W eU nJ\f implies ||^M^||oo= (^W^)(0) and Holder's inequality provides 



/ 1/2 



iAW)\0) 



W{ip) dip 



\-l/2 
I 1/2 



< 




WcL{v)W{ip)dip 



1/2 



W (iff Alp I = / W{LpfALp 



-1/2 



WQL(ip) d(y9 

V-1/2 I Vl/2 

< 1 = (^TycL)'(O). 

Moreover, the estimate in (25) is strict unless there exist a constant c such that 

W{ip) = cW"cl('p) for almost ah e [-i, i], 
whereas the estimate in (26) is strict except for 

1/2 

W{ipfdif = 1. 



(25) 
(26) 

(27) 
(28) 



-1/2 



Now suppose that both (27) and (28) are satisfied. Then we have c = 1, and the norm constraint 
\\W\\< 1 implies W{ip) = = 1^cl(<p) for almost ah tp ^ [-i, i], and hence W = W^cl- □ 

Lemma 3.4. The sequence (*i>n)„ has the complete localisation property on [0, 1] provided that the fol- 
lowing two conditions are satisfied. 

1- VniWcL) ^ I for all n. 

2. $„ converges uniformly and essentially monotonically to on each interval [0, ro] with < ro < 1. 
That means for any ro we have supQ<;^<;^j^ ^nir) — > as n 00, and there exists no(ro) such that 

0<$„,(r)<$„,(r) 

for all 71,2 > Til > no(ro) and all < r < tq. 
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Proof. First we assume additionally that Wn Woo weakly in L^, and suppose for contradiction that 
Woo 7^ WcL- Then, || AW^oo ||oo< 1 thanks to Lemma 3.3, and since AWn AWoo point-wise and 
II AWn ||oo= {AWn)iO) we find some < tq < 1 such that || AWn ||oo< ''o for almost all n. Therefore, 
'Pni {Wn2 ) > 'Pn2 (Wn^ ) holds for all ^2 > Tii and all sufficiently large ni , and since each Wn is a maximisers 
for Vn we also have Vn^iWui) > 'Pn^iWch) > 1- We conclude that P„i(W„2) > 1 and passing to the 
limits 712 ^ oo and rii ~> oo we obtain liminf„^oo ^n(W^oo) > 1- However, || AWoo ||oo< ''o < 1 implies 
lining ooVn (Woo) = 0, the desired contradiction. The result obtained so far implies that Wql is the 
unique accumulation point of a maximising sequence, and this yields the weak convergence to M^cl for 
any maximising sequence. Finally, the strong convergence follows from 1 = ||IVcl||2 =|| Wn II2 for all 
n. □ 

Notice that the two conditions from Lemma 3.4 imply <f>,i(r) for all < ?' < 1, but 00 
as n — > 00, and for this reason it is not clear whether or not the complete localisation (convergence 
of maximisers) implies the convergence of maxima. The simulations from Figure 6, however, provide 
evidence for liminf 7'„(iy„) > 1 = lim VniWch)- 

n — »oo n — ^00 

Our first application concerns the maximiser for homogeneous potentials of large degree. For an 
illustration we refer to the numerical results in Figure 5 and 6. 



Example 3.5. The family of potentials ^q{r) 
property on [0, 1] for q ^ 00. 



;-!-r* with c > and q > 2 has the complete localisation 



Proof. Vq{WcL) = 1 follows from (24) by a direct computation, and for fixed < ro < 1 wc choose qo 
such that l + {qo + 1) Inro < 0. Then we find 9g(^r«) = ir«(l + (g + 1) Inr) < for aU < r < ro 
and q > qo- □ 
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Figure 5: Profile functions W with 7 = i, L = 2, and potentials <i>g as in Remark 3.5. This example 
describes the wave trains for homogeneous potentials in the limit of increasing degree. 
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Figure 6: Speed a and potential energy V{W) versus q for the wave trains from Figure 5. 
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The second candidate for the complete localisation property is related to the limit 7 ^ oo for fixed 
potential For this purpose we consider the rescaled potentials 

<I>,(r) := ^ ^ ^ , (29) 

2 / $(^275) ds 


with corresponding energy functionals V-y, and notice that the two optimisation problems 

V max on Sy and max on S1/2, 

are equivalent due to W £ S-y <^ Wj^flq G 5i/2. 

Example 3.6. For $ as in (22) the rescaled potentials $^ from (29) have the complete localisation 
property on [0, 1] for 7 ^ 00. 

Proof. V^iWch) = 1 holds by construction, and for each < tq < 1 one can find 70 such that d^^.y{r) < 
for all > r > 7'o and 7 > 70- □ 

More generally, the family {^-y)^ can be expected to have the complete localisation property for 7 — > 00 
provided that <I> grows faster than every polynomial. The super-polynomial growth condition is necessary 
as for every homogeneous potential of degree q we have = $1. This refiects the the homogeneous 
scaling 

W-^XW, u'^ ^ Xi-^a^. 

and shows that the wave trains for homogeneous or polynomial potentials do not localise in the limit 
7 00. 



4 Solitons 

This section we study soliton solutions to (12). that means we set 

L = oo, A^A, S^Ur\N. 

Moreover, we assume that the potential energy V is genuinely super-quadratic (see Definition 4.4 below) 
and show that for each 7 > there exists a maximiser of 7^ in 5-,, = U D Af D B^, which is a soliton 
according to Theorem 2.3. 



In what follows we set 



P(7) := sup nW), 



and in order to compare with the harmonic case we introduce 



-Pharm(7) SUp Pharm(VF), (30) 



where 

narm(VK) = / *harm ((^VF) (^)) dip = \(3 \\ AW \\ 



2 



denotes the energy functional corresponding to <&harm('') = ^f^f"^ ■ 

First we show Pharm(7) = Pi by studying the maximizing sequence (!/«)„ C defined by 



t/„(^).| ^-«(^^) M^"' (31) 



for > n. 



16 



Lemma 4.1. We have (i^ > T'harmCC^n) > Pli}- ^ 0{n ^)) for all n, and hence 

Phmm = SUpT'harmU, = SUp Pharm I = Pi- 

Proof. A direct calculation shows i II il2= 7^ ''^cll as 
and 

< {AUM < t/„(n - i) = ^ cos (|(1 - ^)) = 0(n~3/2) ||^| _ „| < 
Moreover, we have 

narm(t/„) >/?7e(£)' 1 COS (|^) ' d^ > /37(l " 0(n-2)) , 



and this implies 7^harm(C/n) — ^ /37 as n ^ oo. Finally, due to Lemma 2.5 we find 

Vh.UW) = i/3 / |(^M^)(^)|' dv5 = /3i II^W^IIB lll^lli 



K 

for all W G B^, and the proof is complete. □ 
Corollary 4.2. We have P{-f) > f3j = Pharm(7) for all 7 > 0. 

Proof. Since the case /3 = is trivial we suppose (3 > 0. Some elementary analysis shows that the 
sequence {Un)n from (31) satisfies 

nUr,)^ J $(e(£)c/„M)d^«i J C/„(^)2(/3 + o(l/V^))d^«/37, 

|ip|<n-i \ip\<n-^ 

where all approximation errors tend to as n — > 00. □ 

As already mentioned in the introduction, solitons are genuinely nonlinear phenomena. In particular, 
the harmonic chain does not allow for solitons, and thus the supremum in (30) cannot be attained. 

Remark 4.3. For /3 > there is no maximiser for T'harm- 

Proof. Suppose for contradiction the existence of a maximizer WVax for T^harm in Bj. Then the Lagrangian 
multiplier rule implies that there exists a (non-negative) Lagrangian multiplier such that 



maxi 



i.e., Wmax is an L^-eigenfunction to with corresponding eigenvalue cr^//?. This is the desired contra- 
diction because the point spectrum of A? is empty due to Lemma 2.5. □ 
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On the weak convergence of unimodal, even and non- negative functions in L^(R) For the 
sake of clarity we proceed with some remarks on the strong compactness of weak convergent sequences 
from Z// n A/'n as this problem becomes relevant in our existence proof for solitons. Strong compactness 
criterions are, in principle, provided by the concentration compactness method from [Lio84], but since 
here we consider only functions from U HAf the arguments simplify a lot. 

Consider a sequence (W„)„ C ^ U Ci H that converges weakly in to some limit Wc^ G S^. 
Passing to a subsequence we can always assume that 7n = ^ || Wn II2— ^ 7oo ^'^^ some 700 with 700 < 7oo ^ 
7, where 700 = ^ II W^oo 111- For 7oo = 7oo the convergence of norms enforces the convergence Wn Woo 
to be strong in L^, and wc are done in this case. 

In the case > joo the convergence cannot be strong as some amount of the 'mass' of the measures 
fin = Wn{ip)'^ d(p disappears when passing to fioo = Woo{^)'^ d(p. However, since all functions are non- 
negative, unimodal, and even, the annihilation of mass is governed by only two elementary processes, 
compare Figure 7. The weakly convergent sequence can form a peak at the origin, and/or a 'fat' tail, 
where 'fat' means 

— L 00 

hm hm / W?,iip)dip+ [ Wl{v)d^>Q, 

— 00 L 

so that some non-negligible amount of the norm is transferred to infinity. 



W„{ip) versus if Wn+mi'f) versus ip 




Figure 7: On the weak convergence in W n A/" n S-y. A weak convergent sequence may fail to converge 
strongly due to the formation of a peak at the origin, and/or a 'fat' tail. 

The first observation is that a peak does not contribute to the potential energy V. In fact, if the 
height of the peak is of order 1/e with e ^ 1, then the norm constraint ^ || ||2< 7 implies that the 
width of the peak is of order e^, and thus the peak disappears after applying the averaging operator 
A. More rigorously. Lemma 2.5 guarantees AWn to converge strongly to AWoo on each compact subset 
of M. Consequently, if the strong convergence fails due to the formation of a peak only, we still have 

V{Wn)^V{Woo). 

The formation of a tail, however, is much more crucial as this in general implies VlWn) ^ 'P(W^oo)- 
Remember the maximising sequence J7„ for T'harm from (31) having the property that all mass of the 
measures Un{f) dip is contained in a fat 'tail' with increasing support and decreasing height, so that the 
weak limit is zero. Even worse, since the spectrum of A^ = /3~^9Pharm is continuous, each maximising 
sequence for Pharm is expected to have this property, and we conclude that the formation of tails is 
directly related to the non-existence of solitons for the harmonic chain. 

Our strategy to prove the existence of solitons for nonlinear potentials is to show that each sequence, 
that maximizes V in Sj, is localised, so a 'fat' tail cannot be formed. To this end we restrict our 
considerations to the case of super- quadratic growth, and derive suitable tightness results under this 
assumptions. 
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4.1 Existence of solitons for genuinely super-quadratic V 

For the remainder of this section we require super-quadratic growth conditions. We start with the 
assumptions concerning the energy functional V as they appear naturaUy in our existence proof. Below 
in §4.2 and §4.3 wc then discuss the corresponding properties of the atomic interaction potential <f>. 

Definition 4.4. The functional V is called super-quadratic on with 7 > if 

V{sW)> s'^nW) for all W e and l<s<j^^. (32) 
Moreover, V is called genuinely super-quadratic on Sj if in addition 

P(7) >^harm(7)=/37- 

Notice that for each 7 > the harmonic functional T^harm is super-quadratic, but not genuinely 
super-quadratic in S-y. 

Remark 4.5. Let V be super- quadratic on S^. Then it is super- quadratic on for a/^ < 7 < 7, and 

we have 

72 ~ 71 ~ 

for all < 7i < 72 < 7- In particular, if V is genuinely super- quadratic for 71 < 7, so it is for every 72 
with 71 < 72 < 7- 

Remark 4.6. Let V be a super- quadratic Oil (S-y . Then, we have 

{dr[W], W) > 2T{W) (33) 

for all VK e 5^. 

Proof. In view of the continuity properties of V and &P it is sufficient to consider the case || W ||2< 7. 
For sufficiently small e (32) implies 

£-^(^'((1 -I- e)W) - V{W)) > £~i((l + e^) - l)r{W) = (2 + e)r{W) 

and the limit £ ^ gives (33). □ 

We now formulate or main technical result concerning the tightness of maximising sequences. Roughly 
spoken, 'fat' tails are not energetically optimal as their contributions to V and T^harm sltc comparable, and 
peaks are not optimal as they do not contribute to the potential energy at all. These naive explanations 
can be stated rigorously as follows. 

Lemma 4.7. For any 5 > Q the set 

S^,s^{W€S^ : ViW)-P^\\W\\l>S} 
is closed under weak convergence. 

Proof. Let (W^n)„ C 5-y be a given sequence such that Wn ^ Woo as n — > 00 weakly in L^, and V{Wn) > 
Pjn + S with 7„ = ^ 1 1 Wn 1 1 2 • Without loss of generality we can assume that ^ 1 1 Wn 1 1 2 — > 7oo some 
with 7oo < 7oo < 7 and 700 = 5 || Woo Hi- It remains to show that V{Woo) > /?7oo + S. For n e N U {00} 
we set 

Wn := Wn\[-m.^+m], Wn := W„ - Wn, (34) 

where m > is some constant to be chosen below, and this definition implies 

\\Wn\\l^\\Wn\\l + \\Wn\\l . (35) 
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Our strategy for this proof is to establish the approximations 

where the approximation error becomes arbitrary small if both m and n are sufficiently large. To show 
this we fix e > 0, and suppose m to be sufficiently large such that 



h \\Woo\\l<e, 



V{W^)-V{Woo) 



(36) 



Such a choice for m exists as VFoo — * Woo strongly in as m — > oo. Since m is finite Lemma 2.5 provides 
AWn AWoo strongly in as n — > oo, and thus we find 



\Woo\\l^\\Wn\\j 



viWoo) - r{Wn) 



< e, 



w^ooM- w-„(v?) d^ 



rn-1 



< £ (37) 



for all sufficiently large n. Moreover, combining (36)i and (37)i with (35) and 7„ shows that 

k II \\l= In - \ II Wn ||^< 7oo - 5 II W^oo ||^ ^Cs < 7oo " 7oo + Ce (38) 

with C independent of n and e. In virtue of (36)3 and (37)3, and since all functions Wn are unimodal, 
non-negative, and even, we also obtain 



< {AWn)W) < iAWn)i^) < iAWn){m-^) < J W„i^)d^ < 2e 

m— 1 

for large n and all (p with \(p\ > m — ^. This provides 

\\AWn\\oo<Ce, J <i>{AW„)dip<Ce, 

I lvl-™|<5 

and exploiting the expansions of $(?■) for small r, compare (14), we find 

V{Wn) < i(/3 + 0(e)) / iAWn)iipf dip < i(/3 + o(e)) || R \\l< f3(-f^ - 700) + Ce (40) 



(39) 



thanks to (38). Finally, due to (34) we have 

r{Wn)^ [ ^AWn)d^+ I <S>iAWn)d^+ I ^AWn)d^ 



\v\<r 



l"Pl>™+2 



<S?{AWn)dip+ / <^>{AWn)d^+ / <P{AWn)dip 



^AWn)dip + riw„), 



|¥'|<m- 2 
< ViWn) + 

|l¥'|-™|<f 

and (36)2 and (37)2 combined with (39) and (40) imply 

ViWn) < V{Woo) + V{Wn) + Ce< V{Woo) + - loo) + Ce. 

Using the assumption we conclude that ViWao) > 7oo + <^ — Ce, and this completes the proof because e 
was chosen arbitrarily. □ 
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As a direct consequence of Lemma 4.7 we find in the genuinely super-quadratic case that each max- 
imising sequence must be localized and contains a strongly convergent subsequence. 

Corollary 4.8. Let V be genuinely super- quadratic on Sj, and suppose that the sequence (W„)„ C S-y is 
a maximising sequence for V on S^, that means 

lim V{W,,) = P(7) =P^ + 5> Pharm(7), 

n — >OQ 

for some (5 > 0. Then, there exists a subsequence, still denoted by (W„)^, and Woo G dS^ such that 
W„ Woo strongly in I?, and hence 'P{Waa) = P{l)- 

Proof. Wc choose the subsequence and Woo G S-^ such that Wn Woo weakly in L^. Thanks to (32) we 
know that V{sWn) > s'^V{Wn) for all s > 1, and this imphes 7,1 = ^ || W„ |||— > 7 as n — > 00, because 
otherwise the sequence (W„)„ could not be maximising. Therefore, with 
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we find V{Wn) -^13^ + 6 and Wn Woo weakly in S-^, and Lemma 4.7 provides 

V{Woo)>Pl^ + 5 

with 7oo = 5 II W^oc lli^i 7- In order to show Woo S dS^ we use again (32) to obtain 

/37 + -5 = P[l) > v[^Woo) > ^ViWoo) > ^{dloo + 5)>p^ + ^S. 

Since 5 > Owe conclude that 700 = 7, that means || Wn Ib^ll Woo II2, and this implies that the convergence 
Wn — > Woo is strong in L^. □ 

The combination of Theorem 2.3 and Corollary 4.8 immediately provides the desired existence result 
for non-negative and unimodal solitons. 

Corollary 4.9. IfV is genuinely super- quadratic on Sy, then there exists a maximiser W of V in S^. 
This maximiser is a non-negative and unimodal soliton with i || WW^— 7. 

Finally, we characterise the soliton speed of maximizers of V , and mention that an analogous result 
holds for wave trains (provided that V is super-quadratic). 

Remark 4.10. The soliton from Corollary 4-9 is super-sonic, that means cr^ > 

Proof Testing the soliton equation (12) with W, and using (33), we find \\W\\l> 2V{W) > 2(3-/ ■ □ 
4.2 Criterions for super-quadratic V 

Definition 4.11. The potential $ is called super- quadratic on the interval [0, \/2^\ if 

$(.sr) > s^$(r) 
holds for all r > and all s > 1 with rs G [0, \/2^ . 



Remark 4.12. // the potential <f> is super- quadratic on the interval [0, \/2'y\, then V is super- quadratic 
on Sj. 

Proof Let W G he fixed, and s be arbitrary with 1 < s < ^27/ || W II2. For all r = {AW){ip) < 
y/\\Wy we have rs < ^27, and hence ^{s{AW){(p)) > s^^{{AW){(p)), for (almost) aU </? e M. Finally, 
integration w.r.t. (p yields the desired result. □ 
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Definition 4.11 implies that $ is super-quadratic on the interval [0, v^^t] if and only if the function 

P^ir) := 2^ (41) 

is non-decreasing on [0, \/2j] (notice that /3$(0) = /3 in the sense of Assumption 2.4). To obtain further 
characterisations of super-quadratic growth we consider the following differential inequalities. 

(CI) $'(r)r - 2$(r) > for all r G [0, 

(C2) $"(r)r - $'(r) > for all r e [0, 

(C3) $"'(r) > for all r e [0, ^2^. 

Remark 4.13. For all j > and all sufficiently smooth potentials $ with $(0) = $'(0) = we have 

{C3) =^ {C2) =^ {CI), 
where ( CI) holds if and only if ^ is super- quadratic on [0, \/2j\ . 

Proof (C3) is equivalent to -^($"(r)r - $'(r)) > 0, and in view of $"(0)0 - $'(0) = we conclude that 
(C3) ^ (C2). Moreover, the imphcation (C2) (CI) can be proven similarly, and (CI) is equivalent to 
■i:P^ir) > 0. □ 

We proceed with a remark concerning the relation between super-quadratic growth and the convexity of 
<f>. In our context, of course, super-quadratic potentials being not convex are forbidden, but solitons can 
still be shown to exist for such potentials, see [FW94] and our comments below. 

Remark 4.14. {C2) implies the convexity of ^ (on the interval [0, y/2^\), but there exists non-convex 
potentials satisfying (CI) 

Proof. Suppose that $ satisfies (C2). Then, the comparison principle for ODEs gives $'(r) > (3r for some 
(3 > 0, and (C2) implies $"(?') > (3 for all r > 0. Now let > be arbitrary, and consider the potential 

$,,(r) = (l + 27r"^ arctan {r]{r - 1))) , 

which is super-quadratic on [0, oo) as the function /3$^ is strictly increasing by construction. A direct 
calculation yields 

(1 + V'^) = vr-i (3 + Stt + 277 - v"^) , 

hence is not convex for large r/. □ 

Another remark concerns the convexity of forces, which become important in the context of atomistic 
Riemann problems in FPU chains, compare [HROSb]. 

Remark 4.15. (C5) implies the convexity of ^' (on the interval [0, ^/Qr/]), but there exists potentials $ 
that satisfy {C2) with non-convex derivative. 

Proof. The first statement is obvious, and towards the second claim we argue as follows. We choose a 
non-negative, but not monotonically increasing function h with /i(G) = 0, and compute $' as solution to 
the ODE h{r) = <I>"(r)r — <i>'(r). Then, each local extremum of h for r > is a turning point of and 
vice versa. □ 
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4.3 Criterions for genuinely super-quadratic V 

In order to complete the existence proof for solitons we must show that for a given super-quadratic 
potential $ the corresponding energy functional V is in fact genuinely super-quadratic on S^. In the 
simplest case there is no harmonic contribution to V at all, and then there exist solitons with arbitrary 
small 7. This holds in particular for all homogenous potentials <I>(r) = cr" with c > and a > 2. 

Remark 4.16. Let $ be super- quadratic on [0, 00) with /3 = $"(0) = and $(r) > for all r > 0. 
Then, V is genuinely super- quadratic on S-y for all ^ > 0. 



The case /3 > is more involved and needs a better understanding of the balance between the harmonic 
and non-harmonic contributions to V. Our strategy in this case is to find a particular function Wq such 
that 7^(1^0) > II III, and this in turn implies the existence of solitons for all 7 > ^ || Wq |||. 

Lemma 4.17. Suppose that $ is super- quadratic on [0, 00) and that the function /3$ from (41) satisfies 
linir^oo /30('') > §/? — |/3*(0). Then, V is genuinely super- quadratic on for all sufficiently large 7. 

Proof. For fixed < £ < 1 and W = \/2jWcl with Wcl as in (23) we find 

1 1 



ViW)^2 J $(v/27s)ds > 27 J /34,(v/27s)sM 


e 1 

>27/3<,(0) J s'-ds + 2jP^isV2^) J .2 a, ^ (^^-^ (1 





and conclude that V{W) > jSj for all sufficiently large 7. □ 
Lemma 4.17 implies the existence of solitons for weakly super-quadratic potentials as for instance 
$(r) = |r^(l + cln(l + r)), $(r) = |r^(l darctan (r)) 
with /3 > 0, c > arbitrary, and d > sufficiently large. 

Next we evaluate the sequence ([/«)„ from (31) and find an existence criterion for solitons that is very 
close to that given in [FW94]. 

Lemma 4.18. Let $ be super- quadratic on [0, 00), and suppose 

$(r) > ^I3r^ +erP 

for all r > with two constants e > and p > 2. Then there exists 70 > such that V is genuinely 
super- quadratic on S-y for all 7 > 7o- Moreover, P ~ or 2 < p < 6 implies 70 = 0. 

Proof. Let 7 > be arbitrary, and consider the sequence {Un)yi from (31). Then, 

V{Un) - nao„(C/„) > eVnxiUn) 



ith 



\ip\<n ^ 



2 

J_ 

'2n 



> 



^{w.)rC-^Y''n / (cos(fv.))M^>cn-^/^. 



1 
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We conclude 



r{Un) > 7'harm(t/„) + CjP/^n'-P/^ > 

for some positive constant c > 0, and according to Lemma 4.1 there exists a constant c > such that 

ViUn) > c-i^l'^n^-Pl'^ + /37(1 - 

Finally, for /3 = 0, or 7 sufficiently large, or— f+p/2<2 and n large, we find 'P{Un) > Pj- □ 

As an application of Lemma 4.18 we find the following existence result for unimodal solitons with 
non-negative W: Let $(r) be analytic with non- negative coefficients, i.e., $(r) = ^(3r^ + X^i^a with 
Ki > for all i > 3. Then $ is super-quadratic on [0, 00) and genuinely super-quadratic for large 7. 
Moreover, if at least one of the coefficients K3, K4, and K5 is positive, then <I> is genuinely super-quadratic 
for aU 7 > 0. 

Moreover, since the traveling wave equation (12) is invariant under the reflection symmetry 

W-^-W, $(r)->$(-r), 

we find also existence results for solitons with non-positive W. For instance, the Toda potential $Toda(?') = 
e^'' + r — 1 is not super-quadratic for r > but the reflected potential $Toda(?') = e*" — r — 1 has solitons 
with arbitrary small 7, compare Lemma 4.18. From this we infer that the Toda chain allows for solitons 
with < 0. 

Finally, we summarise some other super-quadratic growth conditions for <f> under which the existence 
of solitons was proved by other authors. 

1. <i>'(r)r > 2$(r): Friesecke and Wattis [FW94] prove the existence of super-sonic solitons with 
prescribed potential energy V > Pq above some critical energy Pq > 0. Moreover, $(r) = i/3r^ + 
£:rP(l -I- o(r)) with e > and p as in Lemma 4.18 implies Pq = 0. 

2. <I)(r) = -f ^ni{r) and $ni(?')r > a^niir) for aU r > and some a > 2: Smets and Willem 
[SW97] establish the existence of solitons with prescribed super-sonic speed cr^ > /? = $"(0) 

3. <I>(ji(r)r > a^niir) for all r > and some a > 2, or $"i(?')r > Q;$J^[(r) for all r > and some 
a > 1: Pankov and Pfliiger [PPOO] prove under both assumptions the existence of a family of 
solitons parametrised by cr^ > /? = $"(0). 

4. <I>'(r) > for aU r > and liminfr^oo r^°'{^'{r)r — a$(r)) > and some a > 2: Schwetlick and 
Zimmer [SZ07] show that for each super-sonic speed > cTcrit > = <I>"(0) there exists a soliton. 

All these existence results imply that the soliton profile W belongs to Af, but since they do not require 
the convexity of they do not provide W £U. 

4.4 Solitons as limits of wave trains 

It is very natural to ask whether or not wave trains converge to solitons when the periodicity length L 
tends to 00. In this section we establish such a convergence result for unimodal and non- negative wave 
trains. To this end we allow for arbitrary values of L G (0, 00], and write Al for the operator A acting 
on L^([— L, L]). Consequently, we introduce 

SL,-y ■.= ur\Ur\BL.-y 

where Bl.-j denotes the ball of radius y/2rf in L^([— i, L]), and consider 

L 

Pl-^ sup [ ^{{ALW){^))d^. 

W&Sl J 

-L 
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Moreover, we define an embedding operator El : Sl,^ Soc.-y by 

Inspired by the notion of r-convcrgence we show that the energy of each periodic profile can be 
approximated by localised profiles, and prove that each localised profile can be recovered by periodic 
profiles. 



Lemma 4.19. For each 7 > and L < 00 there exists a constant C^^^ of order o{yJ ^ / L) such that 

Voo{ElWl) + Cl,^ > Vl{Wl) > Voo{ElWl) (42) 

holds for all Wl G SL,-f Moreover, for any Woo G 500,7 there exists a family of functions {Wl)l^oo 
such that 

Vl{Wl) Voo{Woo)- (43) 

Proof. First let Wl G Sl.^i be fixed, and notice that 

{AooElWl){.v) = {ElAlWl){v), \\v\~L\>\, (44) 
hold by construction. For 1 — L| < i we have < Wl{,v) < Wl{L - 1) due to Wl 6 UnAf, and hence 

< {AooElWl){^) < (ElAlWlKv^) < J WL^d^ = Wl{L-1). 

Moreover, from Wl eUnU we infer {AooElWl){(p) < {ELALWL){f) and 

L-l 



27 > J WLi^f > 2{L - 1){Wl{L - l)y 



-L+l 

and therefore 



< iAooELWL)iip) < {ELALWL)iip) < e, \\^\ - L| < i. (45) 



< / $ 



with e = \J^I(L — 1). The estimate (42) now follows from (44) and (45) via 
^<Vl{Wl)-Voo(ElWl) 

\ElAlWl){'^)) - ^({AooElWl){'^)) d(^ 
< y" ^{{ElAlWl)W)) - '^({AooElWl){'^)) d(/? < /(e), 

||¥>|-l|<1/2 

with /(e) = 2c sup $'(r) = o(e). Now let W^o G iSoo,7 be fixed, and define Wl & Sl,^ by 

0<r<e 

{Wl){^)^Woo{^) for 1^1 <L. 

Then, in general we have Wooif) 7^ {ELWL){f) = for > L but always ElWl Woo strongly in 
L2(R) as L ^ c». This implies VooIelWl) ~ 'Poo(M^oo) 0, and thanks to (42) we find (43). □ 
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Lemma 4.19 now provides the convergence of suprema. 
Corollary 4.20. We have Pl{i) -'^^^ Poo{l)- 

Proof. For given L < oo let W l be a maximiser of V on Sl.i- Then (42) implies 

Pdl) < VooiEiML) +oiVl/L) < Poo (7) + 0(^77^), 
and hence limsup^_,o^ Pi(7) < Poo (7)- Moreover, in view of (43) we have 

Poo(Woo) <liminfPi(7) 

L — ^00 

for all M^oo, and this shows Poo (7) < liminf^^oo Pl{i)- D 

As a further consequence we find that solitons can be constructed as limits of wave trains. More 
precisely, Corollary 4.20 combined with Corollary 4.8 provides the following convergence result for max- 
imisers. 

Corollary 4.21. Let V he genuinely super- quadratic on S^, and for each L < 00 let Wl be a maximiser 
of V in Sl,^- Then, for any sequence (iji)„ with L„ — > 00 there exist a subsequence, still denoted by L„, 
and a maximiser Woo G 500,7, such that El^Wl„ Woo strongly in L^(E). 
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